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We describe a broad class of bounded non-periodic potentials in one-dimensional stationary quan¬ 
tum mechanics having the same spectral properties as periodic potentials. The spectrum of the 
corresponding Schrodinger operator consists of a finite or infinite number of allowed bands sepa¬ 
rated by gaps. In this letter we consider the simplest class of potentials, whose spectra consist of an 
interval on the negative semiaxis and the entire positive axis. The potentials are reflectionless, and 
a particle with positive energy moves freely in both directions. The potential is constructed as a 
limit of Bargmann potentials and is determined by a Riemann-Hilbert problem, which is equivalent 
to a pair of singular integral equations that can be efficiently solved using numerical techniques. 
PACS numbers: 02.30.Ik, 03.65.Nk 


Introduction. We consider the Schrodinger equation on 
the real line: 

— 'll)" + u{x)ij) = Elf for —oo < X < oo, (1) 

where u{x) is a bounded potential. The spectrum of u{x) 
is the set of values of E for which there exists a bounded 
solution ijjlXjE) of (P). 

Traditionally, the goal of spectral theory in quantum 
mechanics has been to describe the possible spectra of a 
given class of potentials of ©• For an arbitrary bounded 
potential, the answer is not known in general. Two cases 
have been extensively studied and are well understood: 
the rapidly decaying case and the periodic case. The 
spectrum of a rapidly decaying potential consists of the 
positive real axis and a finite discrete set on the negative 
real axis. For a periodic potential, the spectrum is purely 
continuous and consists of an infinite collection of allowed 
bands separated by spectral gaps. Potentials having only 
finitely many spectral gaps form a dense subset among 
all periodic potentials. These potentials are described 
by algebro-geometric methods and are expressed by the 
Matveev-Its formula in terms of the Riemann 0-functions 
of hyperelliptic curves. The Matveev-Its formula also 
describes a class of quasi-periodic potentials having a 
finite-gap structure. However, the spectrum of a generic 
quasi-periodic potential has a complicated Cantor set-like 
structure, and the corresponding system displays Ander¬ 
son localization. 

In this letter, we formulate the opposite problem of 
describing the set of potentials of the one-dimensional 
Schrodinger operator that have a given spectrum. Specif¬ 
ically, we pose the question of describing the class of 
bounded potentials of dm having the same spectrum as a 
finite-gap algebro-geometric potential, but which are nei¬ 
ther periodic nor quasi-periodic. Physically such a poten¬ 
tial describes a non-periodic medium which nevertheless 
admits free wave propagation. Solving this problem is 


especially relevant given the recent explosion of research 
into metamaterias that have predetermined optical prop¬ 
erties or a given electronic band structure. 

In this letter we focus our study on one-gap potentials, 
with a spectrum consisting of the positive real axis and 
one allowed band [—k 2 ,—kf] on the negative real axis, 
where 0 < ki < k 2 . A periodic one-gap potential with 
such a spectrum is determined by the formula 

u{x) = 2p(x + iu)' — xq) + 63 , (la) 

where p(z) is the elliptic Weierstrass function with peri¬ 
ods 2uj and 2iuj', and 

kj = ei - 63, fci = 62 - 63, 61-1-62+63=0, (2) 

where ei > 62 > 63 are the values of p{z) at the half¬ 
periods of the lattice. The spectrum of the correspond¬ 
ing operator is doubly degenerate and reflectionless. We 
show that a generic one-gap reflectionless potential is de¬ 
termined by two positive continuous functions defined on 
the allowed band. 

Bargmann potentials via dressing method. We describe 
an efficient method for constructing non-periodic one-gap 
potentials by taking the closure of the set of reflectionless 
Bargmann potentials [I| . This construction was formally 
carried out in the works of Marchenko et al i,i and Q, 
but their results are not effective. We support our results 
with numerical simulations. 

Since the discovery of the 1ST in 1967 the 
Schrodinger operator has been used for exact integration 
of the Kortweg-de Vries (KdV) equation. The reflection¬ 
less Bargmann potentials were reinterpreted as the time 
frames of the fundamental A^-soliton solutions of KdV. 
For this reason, and also because the dressing method is 
local in time, we retain this terminology. 

The most direct and elegant method for constructing 
Bargmann potentials, known as the dressing method, was 
described in [^. The dressing method is a considerably 
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more flexible tool than the traditional Inverse Scattering 
Method (ISM). 

Following the cited paper, we consider the following 
i9-problem on the complex fc-plane: 


The corresponding potential 


N 




(9) 




( 3 ) 


where T{k) is a compactly supported distribution called 
the dressing function of the 5-problem. The 5-problem 
is normalized by the condition % —>■ 1 as |A:| —>■ oo. The 
solution satisfies the following integral equation: 


TT i/ k + q 


where we regularize the integral in the following way: 


^ = lim . E „ ■^l=TTS{k), 


k e^-o dkk 


is a reflectionless Bargmann potential. The function 
Ip = satisfies the Schrodinger equation ([1]). The 

corresponding potential is rapidly vanishing and has a 
finite discrete spectrum {—,..., —k^}, and ipn{x) = 
Xn(x)e'^”^ are the corresponding eigenfunctions. 

The basic relation is equivalent to the following 
linear system on the eigenfunctions: 


tpnix) -i-T„y^ 




n—1 


“ 1 “ 


— TnC 


( 10 ) 


iiTn/ Kn > 0 the determinant of this system is positive 


A = det 


> 0 , 


Kn T 

SO system cni) has a unique solution, and moreover 


( 11 ) 


where 6{k) is the two-dimensional J-function. Suppose 
that the dressing function T{k) has the property that 
equation ([I]) has a unique solution. Then 


x{x, fc) —^ 1 -I- -I-... as |fc| —>■ oo (5) 

k 


and the function x is a solution of the equation: 


Xxx-‘2ikxx-u{x)x = 0, u{x) = 2—xoix). ( 6 ) 

ax 


To construct Bargmann potentials, we consider the dress¬ 
ing function 


r(fc) = 7r^T„d(/c - iKn), 


( 7 ) 


where Tn and are nonzero real numbers satisfying 
\txn\ G [^ 17 ^ 2 ] and TnjKn > 0 for all n. Then x is a 
rational function: 

xix, k) = 1 + i X’^(^) where Xn{x) are real. ( 8 ) 


u{x) = - —^\nA 
ax^ 


( 12 ) 


Closure of Bargmann potentials. To describe the clo¬ 
sure of the set of A^-soliton solutions, we replace the fi¬ 
nite dressing function © with one supported on two cuts 
[ikiAk^] and [—ifc 2 , —iki] on the imaginary axis: 


pk2 pK,2 

T(fc) = TT / Ri{K)S{k—iK) dK+TT / R2{K)S{k+iK) dn. 
J ki J ki 

(13) 

Here Ri and R 2 are functions on [fci, k 2 ]. Formula dl]) 
transforms into the following spectral representation: 


r-fc2 


= + ( 14 ) 

ifci k-ip Jki k + ip 

where f{x,p) and g{x,p) are real-valued functions. Then 


rk2 

Xo= [f{x,p) + g{x,p)]dp and u(x) = 2 
Jki 


dxo 

dx 


( 15 ) 

The function xix,k) is analytic away from the cuts. The 
details of the transition from the discrete to continuous 
case are described in Q. The algebraic system (ITOl) trans¬ 
forms into a system of two integral equations: 


f{x,k)+Ri{k)e-^'^^ 
g{x, k) + R 2 {k)e^^'^ 


pk2 

Jki 


i 


/(^,g) 

/Cl k-q 



r f’”’ 

Jkx k-q 

= i?i(fc)e-2'=“ 

(16) 

dq+ 1 
J 

C, k + q 

= -R2{k)e^^^ 

(17) 
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In 0 we show that if i?i and R 2 are continuous, positive 
and a-H61der for a positive a, then system msD-dni) has 
a unique solution. Moreover, 

k) = f{x, and ip{x, k) = g{x, 

( 18 ) 

are eigenfunctions of the Schrodinger operator © cor¬ 
responding to E = kF'. These eigenfunctions form an 
orthogonal collection: 


where S{k) is an odd function (for example, for KdV 
S[k) = 8 k^t). A periodic potential remains periodic 
at any moment of time, so applying transformation (12411 
does not change the periodicity properties of the poten¬ 
tial. We illustrate this with a numerical simulation. 

We now construct a dressing function of the form (fT3l) 
that determines a periodic potential. We put Xq = ui 
in dial) and map the /c-plane to the period parallelogram 
as follows: 



1 

II 

(19) 


pOO 

/ 1 p{x, K)'lp{x, k) dK = R 26 {k — k), 

J —00 

(20) 


0 

II 

(21) 

On any interval where one of the functions Ri, R 2 is 
positive and the other vanishes, the spectrum is simple. 
If Ri and i ?2 both vanish, the corresponding segment is 


a spectral gap, so the finite-gap case is included in our 
description. If Ri = i? 2 , then u{x) = u(—x) and the 
potential is symmetric. 

We also show in Q that system (HHi-dni) is equivalent 
to the following Riemann-Hilbert problem on x- 


k^ = 63 — p{z) z{k) —>■ —— as |A:| —>■ 00 . (25) 


Then ([T]) becomes the Lame equation 

p" — [2p{x — u}' — iu}') + p{z)] ip = 0 (26) 

which has following solution: 

, , a(x — uj — iu + z)<t{u} + iuj') , x ^ 

^) = — —T—-:t exp {-C{z)x). 


a{x — u — iu)a{u + iu' — z) 
Now define the function: 

^{x,k) = ^-^if{x, z{k))e''^^. 


(27) 

(28) 


X'^ix,ik)-x ix,ik) = 

= iTrRi{k)e~‘^^^ [x^(^) —ik) + x~i.x, —ik)] , 


( 22 ) 


This function satisfies the equation: 

— 2ik^' — u{x)^ = 0, for ^ t 1 |fc| —>■ 00 (29) 


X+(x,-ifc)-x {x,-ik) = 


(23) 


- *7ri?2(A:)e^'"^ [x^{x, ik) + x {x, ik)] , 

where x^ are the boundary values of x along the cuts: 

X^ix, k) = lim xix, k±ie), k G [—ik 2 , — 1 ^ 2 ]. 
£^0 


The function ^ satisfies the RH problem (j2^ - (f23|) . where: 
Riiq) = -h{q), i? 2 (g) = (30) 


, . N {k2-q){q + ki) , ^ ^ v. 

'■(«) = V (,-t,)(,+t,) ' ^^ 


A fundamental difference between the ISM and the 
dressing method is that in the former, the scattering data 
can be uniquely restored from the potential, while in the 
latter, the same potential can be constructed using a va¬ 
riety of different dressings. For example, in t © we can 
change the signs of one of the T„, and the corresponding 
without changing u(x) (see Q for details), so each 
A^-soliton solution can be constructed using 2^ different 
dressing functions of the form © . Hence we cannot ex¬ 
pect a given bounded potential to be determined by a 
unique choice of i?i and i? 2 - Below we describe a class 
of dressings leading to one-gap periodic potentials. 

Periodic potentials. Suppose that a potential u{x) 
evolves according to some equation of the KdV hierar¬ 
chy (see Q). The dressing functions evolve as follows: 


Now S{q) = lnh{q) is an odd function, so we can replace 
h{q) with I and assume that 

R3iq)=R2iq) = -. (32) 

TT 

The dressings dSUll - dHT]) and (15^ both define a periodic 
potential. 

Numerical solution. We solve equations dUD-dni) nu¬ 
merically for ki = 2 and ^2 = 4. It is convenient to 
replace ^(x, k) and ip{x, k) with the following functions: 

P(x,k) = \/l — k‘^4){x, k), 

Q{x, k) = \/l — ki^tp^x, k), 


R 2 {k) ^ R 2 {k)e (24) which then satisfy the following integral equations: 


Ri{k) Ri{k)e‘ 
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FIG. 1. A potential U{x) that appears as a result of dressing with: (a) Ri — R 2 = I/tt (periodic); (b) Ri = I/tt and R 2 = 0 
decaying for a; > 0; (c)i?i = ^ x 10“® and R 2 ~ ^ x 10“®. 


P{x, k) + ri{x, k) 
Q{x, k) + r 2 {x, k) 


' P{x,q)e-^^dq 

Q{x,q)e^^dq ' 

Jki {k + q)^/l-q'^ J 

ki {k - q)^/l - q'^_ 

P{x,q)e-‘>=^ dq 

Q{x,q)e^^dq ' 

Jki (fc-g)v^l-g2 J 

ki ik + q)y'l-q'^_ 


ri(a;, fc), 


-r2{x,k), 


(33) 

(34) 


where rj{x,k) = \/l — kP'Rj{k)e^ lykx £gj. j _ ^ 2_ 
Discretized at Chebyshev nodes qk = cos 
k = 1,2,..., M the integrals are evaluated via Gauss- 
Chebyshev quadrature that is exact for polynomials of 
degree less than 2M — 1. Note that each equation of 
the system contains a Cauchy principal value integral 
denoted by f, and that integration in the vicinity the of 
singularity a,t q = k requires a shift from real axis. 

The spatial variable x appears as a parameter in (j33ll - 
(|M1) and the cc-dependence of ri and r 2 becomes a major 
obstacle, since the condition number of the discretized 
system is exponential in x and requires usage of multi- 
precision arithmetics. 

A mesh for the spatial parameter cc is a Chebyshev grid 
with M nodes; a high-order polynomial interpolation by 
means of Lagrange interpolation is used for intermediate 
points. In a typical simulation an interpolating polyno¬ 
mial of degree 200 suffices to have an accurate approxi¬ 
mation for \x\ < 10 (See figures [T|). 

Conclusion. The procedure of closing the set of 
Bargmann potentials described above can be general¬ 
ized to a wide class of linear operators, such as the 
Dirac operator. This is especially important for vari¬ 


ous nonlinear integrable systems, such as the Nonlinear 
Schbdinger Equation and the Kadomtsev-Petviashvili 
equation, where these linear operators play a fundamen¬ 
tal role and which are solved using the ISM. 

To apply the theory that has been developed for these 
systems to real physical problems, it is necessary to de¬ 
velop a statistical theory of integrable systems with an 
infinite number of degrees of freedom. This theory of in¬ 
tegrable turbulence is still very much in its infancy, the 
first stages have been suggested in Q . The technique de¬ 
scribed in this letter shows an approach to constructing 
strongly nonlinear statistically homogeneous solutions 
to integrable systems such as KdV and the nonlinear 
Schrodinger equation. In fact, integrable turbulence is 
a common physical phenomenon. This turbulence takes 
place in the coastal areas of seas, and describes effects oc¬ 
curring in optical fibers. Thus the natural extension and 
our next step will be the study of random non-periodic 
potentials with a continuous spectrum, with this work 
being at the core of the theory of integrable turbulence. 
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